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Abstract 

An elastic anomaly, observed in the heavy fermi liquid state of Ce alloys (for 
example, CeCug and CeTe), is analyzed by using the infinite-f/ Anderson 
lattice model. The four atomic energy levels are assumed for f-electrons. 
Two of them are mutually degenerate. A small crystalline splitting 2A is 
assumed between two energy levels. The fourfold degenerate conduction 
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bands are also considered in the model. We solve the model using the mean 
field approximation to slave bosons, changing the Fermi energy in order to 
keep the total electron number constant. The nonzero value of the mean 
field of the slave bosons persists over the temperatures much higher than the 
Kondo temperature. This is the effect of the constant electron number. Next, 
the linear susceptibility with respect to A is calculated in order to obtain the 
renomalized elastic constant. The resulting temperature dependence of the 
constant shows the downward dip. We point out the relation of our finding 
with the experimental data. 
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1. Introduction 



Generally speaking, the atomic f-level of heavy fermion Ce compounds is 
split by a crystalline field into multi sublevels [1]. The original f-level of 
the degeneracy 14 is first split into the j = 5/2 and 7/2 levels by the l-s 
interaction. The j — 7/2 level is about lO^K higher than the j — 5/2 level, 
so the main contribution comes from the j = 5/2 level, and the j = 7/2 level 
can be neglected. Secondly, the j — 5/2 level splits into several multiplets 
in the presence of the crystalline field. For example, in the cubic field, the 
j — 5/2 level is composed of two sublevels of degeneracy two and four, and in 
the tetragonal field it sphts into three Kramers doublets. When the ground 
level is doubly degenerate, the lowest f-level is nearly half filled because the 
number of the 4f-electron is close to 0.9. The splitting width is normally 
10^~^K, so the temperature dependence of the physical quantities refiects 
the effect of the thermal excitations among the sublevels. 

The effect of the crystalline field on the single site Kondo system has been 
one of recent topics in the theory of dilute magnetic alloys. The Anderson 
model or the Coqblin-Schrieffer model, with the inclusion of the crystalline 
field and the anisotropy of the mixing interaction, is solved by the use of the 
Bethe ansatz technique [2] and the self-consistent renomalization [3]. The 
result shows the variation of the Kondo temperature Tk with the temper- 
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ature. It comes from the change of the impurity scattering channel. For 
example, we imagine the case where there are three Kramers doublets with 
the splitting Ai amd A2. In the temperature range T < Ai, the three levels 
act independently in the scattering. But, in Ai < T < A2, the lower two 
levels behave as if they are one level of the degeneracy four, because thermal 
fluctuation has a larger width than Ai. In the same way, in T > A2, all 
the three levels look as a level sixfold degenerate. It has been discussed that 
the above features are exhibited in temperature dependence of the magnetic 
susceptibility. 

Experimentally, the heavy fermion systems show elastic anomalies related 
with valence instabilities in low temperatures. The most striking effect is the 
softening of the elastic constants below the Kondo temperature which are 
observed prominently in CeAla [4], CeCue [5], and CeRu2Si2 [6]. The tem- 
perature dependences of the elastic constants have been theoretically fitted 
quantitatively well [7] by using the electron phonon coupling derived from 
the ansatz that the mixing strength depends on the volume of the crystal. 
The same ansatz was highly successful in the theoretical description [8,9] of 
the Kondo volume collapse transition between a- and 7-Ce. These elastic 
anomalies are related with the overall instability of the valence mainly due 
to the change of the mixing interaction. They are observed in the wide tem- 



-4- 



perature range < T < Tk and rather insensitive to the detailed structures 
of the electronic band structures. 

However, when we look at detailed temperature dependences of the elastic 
constants, several compounds show anomalies which have been interpreted 
as the crystalline field effect. For example, the elastic constant C33 of CeCue 
[10] has a dip at about lOK. This might be due to the sphtting larger than 
the Kondo temperature 4K. The constant (Cn — Ci2)/2 of CeTe [11] shows 
the apparent dip at about 15K. The ground state levels of Ce ions split into 
the Fy Kramers doublet and the Fg quartet states. The Fy states are the 
ground states. There is the splitting 30K between F7 and Fg states. This is 
the origin of the dip. 

The main purpose of this paper is to present a microscopic calculation in 
order to look at how the degeneracy structure and crystalline field appear in 
the elastic properties. The formalism is independent of the real band struc- 
tures and degeneracy structures. Therefore, the results should be general 
enough for heavy fermion systems. We simply assume four quantum num- 
bers for f-electrons. Two of them have the same atomic energy level. Thus, 
two different atomic levels are assumed. The difference of them, smaller than 
Tk, means the width of the crystalline field splitting. And, four conduction 
bands are assumed. The total electronic system is described by the infinite- 
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U Anderson lattice hamiltonian by the slave boson method. The model is 
solved by the mean field approximation. The formalism is explained in §2. 

Firstly, we show mean field solutions and see that the mean field value of 
slave bosons does not vanish even at T > 7k- This interesting property has 
been discussed previously [12,13]. It is owing to the fact that the Fermi level 
is varied in order to keep the total electron number constant. This is shown 
in §3. 

Secondly, we calculate the linear susceptibility with respect to the crys- 
talline field splitting. At lower temperatures, the susceptibility has a struc- 
ture related with the crystalline field effect. The high temperature suscep- 
tibility does not depend on the splitting width due to the large thermal 
excitation. There is an upward dip in the temperature dependence. The 
dip is originated from the high degeneracy, and its position moves with the 
splitting width. 

Next, we assume an empirical relation between an elastic constant and 
the susceptibility. The form of the relation is assumed by taking account of 
the quadrapoler response theory [14]. As the coupling constant between the 
lattice and the crystalline field is unknown, we should treat it as a kind of a 
parameter. We show temperature dependences of elastic constant for several 
choices of the coupling. We will discuss relevant parameters for the elastic 
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anomaly, i.e., the downward dip which is observed in the constant (Cn — 
Ci2)/2 of CeTe [11]. The susceptibihty and elastic constant are reported and 
discussed in §4. 

Finally, the paper is closed with several remarks in §5. 

§2. Formalism 

We consider the infinite-[/ Anderson lattice model in the slave boson method. 
The model has the following form: 

i 1=1,2 1=3,4 

+ E 

k,l=l-4 

+ V y: ifiiCi,ibi+bi4,ifi,i} 

jJ=l-4 

+ EM E fliAi + blh-i), 

i 1=1-4 

where fi^i is an annihilation operator of the f-electron of the l-th orbital at the 
i-th site, ^ is an operator of the conduction electron with the wave number 
k, and bi is an operator of the slave boson which indicates the unoccupied 
state at the f-orbital. The atomic energy of the first and second orbitals of 
f-electrons is Ei — A, and that of the third and fourth orbitals is Ei + A. 
The crystalline field splitting 2A is assumed between atomic energies of the 
two groups of f-electrons. For the conduction electrons, the same quantum 
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number is assumed as that of the f-electrons. We use the square density of 
states, p = 1/ND, which extends over the energy region, —D < < {N — 
l)D, where = 4 is the total number of quantum states. This assumes that 
the combination NpV^, which appears in the 1/N expansion, is independent 
of N. Therefore, the mean field theory becomes exact as N ^ oo. The third 
term in the hamiltonian is the mixing interaction between f- and c-electrons, 
V being the interaction strength. The last term limits the maximum number 
of f-electrons per site up to unity. This could be realized by the constraint 
J2i=i-4 flifi,i + blbi — 1 with the Langrange multiplier field Aj. 

This model is treated within the mean field approximation: {bi) = r, 
(blbi) — r^, and Aj = A (a site independent real value). These mean field 
parameters are determined by solving the following coupled equations: 
(1) the constraint condition. 



2D J 



dE^ f(E - n) (2) 

1 r V'^ 

+ — dE^ f(E-pL)+r^^l, 

2D J {Ei + A-Ey-'^ ^' 

(2) the self-consistency condition for r, 

r2 



If 

+ — dE . /(£; - /x) + A = 0, 

2D J E-Ei-A^ ' 
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and (3) the conservation condition of electron number riei, 



(4) 



(^f - A - E) 



(^f + A - E) 



■^]f{E - IJ)= ^el, 



where f{x) = l/[exp(a;/T) + l] is the Fermi distribution function, E^ = E{ + X 
is the effective f-level, and V — rV is the effective mixing interaction. The 
integrations are performed over aU the energy region of the bands. The three 
equations are solved numerically for the variables, r. A, and the Fermi level 
/I. In addition, the values at T = can be obtained analytically. 



Equations (2), (3), and (4) are solved numerically for the parameters D = 
5 X lO^K, V = 7500K, Ef = -lO^K, and n^i = 1.9 as the typical values. The 
splitting parameter A is increased from but the magnitude is taken below 
the Kondo temperature — Ef — fj,. 

Figure 1 shows the temperature dependences of parameters. Figures 1 
(a), (b), and (c) show the variations of Ef, Tk, and the number of f-electrons 
per site nf, respectively. The data for A = 0, 15, and 30K are shown by the 
dashed, thin, and thick lines, respectively. This convention applies to all the 
figures in this paper. As the temperature increases, the order parameter r 
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decreases, so that Uf = 1 — increases. The quantity r does not vanish even 
though the temperature is much higher than Tk at T = 0. This is the effect 
of the change of the Fermi level /j, to keep the total electron number constant. 
This effect has been reported previously [12,13]. According to the increase of 
rif, Ef decreases, which means the reduced itinerancy of f-electrons. At high 
temperatures, T ^ A, the three curves becomes almost identical, owing to 
the large temperature excitation across the split 2A. At low temperatures, 
the excitation energy is limited by the smaller distance from the Fermi level 
to the gap of the bands I — 1,2. This results in the increased value of rif 
when the crystalline field is switched on. 

We can derive the analytical expressions of parameters at T = 0. Small 
crystalline field splitting is assumed for the derivation. The results are 



where Tk(0) = D exp[—D{ii — Ei)/V^] is the Kondo temperature for A = 0. 
These expressions accord with the numerical results that both Tk and Uf 
increase as A increases. 




(5) 



and 



V^TKiAY 



(6) 
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§4. Elastic anomaly in low temperatures 

We shall discuss the change of elastic property of heavy fermions due to the 
crystalline field splitting in the low temperature below Tk- Generally, an 
elastic constant c is related with the linear susceptibility with respect to A, 
as shown by the formula [14]: 



where Cq is the elastic constant of the system where there is not interactions 
between the lattice and the electronic system, and g is the coupling constant 
between A and the strain field. We assume that cq is independent of the 
temperature. The value of g is unknown experimentally as well as theoreti- 
cally. In order to discuss the crystalline field effect on c, we treat the factor g 
as a kind of fitting parameters. The quantity xa is calculated as the second 
order derivative of the mean field free energy: 



where the A dependences of the band edges are neglected in the derivatives 
because their effect is exponentially small. 



Xa = 



9A2 



(8) 
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Figure 2 (a) displays the temperature dependence of xa- The inverse of 
Xa is shown in Fig. 2(b). There is a peak at the low temperature. When 
A = 0, Xa is identical with the magnetic susceptibility. The peak of the 
curve for A = has been found in the previous paper [12], and is the effect 
of the fourfold degeneracy. When A is finite, the peak moves to the lower 
temperature and the peak height becomes taller. This is the crystalline field 
effect. The curve at high temperatures is less affected by A. The inverse of 
the susceptibility in high temperatures is almost linear against T, showing 
the Pauli behavior. The increase of xa in the temperatures below the peak 
is well explained by the analytical expression: 
^ Tg(0) + 2A^ 

at T = 0. 

We show the temperature dependences of c/cq, assuming several values 
for g. We plot two sets of curves for g = IK and 1.5K, in Figs. 3(a) and 
(b), respectively. The elastic constant decreases from much higher to lower 
temperatures than Tk- This is the effect of the valence fluctuation. There is 
a downward dip around 15-20 K, which is the result of the combination of 
the fourfold orbital degeneracy and the crystalline field splitting. The larger 
A gives rise to the larger dip. 
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We shall look at relations with actual compounds. The elastic constant 
C33 of CeCue [10] has a dip at about lOK. This might be due to A larger than 
the Kondo temperature 4K. The constant (Cn — Ci2)/2 of CeTe [11] shows 
the apparent dip at about 15K. The ground state levels of Ce ions in CeTe 
split into the Kramers doublet and the Fg quartet states. The F7 states 
are the ground states. There is the sphtting 2A — 30K between F7 and Fg 
states. Fig. 1 of ref. [11] shows that the decrease of (Cn — Ci2)/2 constant 
from that of the ideal system without the electron-crystalline-field coupling 
is about 4% of the magnitude. This behavior is simulated well by the thin 
curve in Fig. 3(b). The value g = 1.5K gives us a useful information on the 
couphng strength of heavy fermions against the crystaUine field in CeTe. 



§5. Concluding remarks 

We have solved the mean field equations of the Anderson lattice model with 
crystalline field splitting 2 A, and calculated the linear susceptibility with 
respect to A. Next, we have simulated the temperature dependence of the 
elastic constant which is derived by the coupling of the electronic system to 
the crystalline field. We have discussed the relation with experiments. 

The downward dip of c exists even if there is no crystalline field split- 
ting, i.e., A = 0. The finite A results in the motion of the dip to lower 



temperatures. The magnitude of the dip becomes larger at the same time. 
The temperature where the dip is present would be finally determined by 
the combination of the high degeneracy of the 4f orbitals and the splitting 
width. 

One of the interesting future problems would be the magnetic field effects. 
Further removal of the degeneracy by the magnetic field will give rise to 
more structures in the elastic constants. The theory of the effects would be 
a helpful information for experimental studies. 
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Figure Captions 

Fig. 1. Temperature dependences of the mean field solution: (a) E{, (b) Tk, 
and (c) rif. Parameters are D = 5 x lO^K, V = 7500K, E{ = -lO^K, and 
Hei = 1.9. The dashed, thin, and heavy lines are for A = 0, 15, and 30K, 
respectively. 

Fig. 2. Temperature dependences of the hnear susceptibihty: (a) xa and (b) 
1/Xa- The parameters and notations are the same as in Fig. 1. 

Fig. 3. Temperature dependences of the elastic constant c/cq for (a) g = l.OK 
and (b) g — 1.5K. The parameters and notations are the same as in Fig. 1. 
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